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Abstract 

So far the nuclear PDFs (nPDFs) in the global DGLAP fits have been taken to be spatially independent. In 
this work [I J, using the A-dependence of the globally fitted sets EPS09 and EKS98, we have determined the spatial 
dependence of the nPDFs in terms of powers of the nuclear thickness functions. New spatially dependent nPDF sets 
EPS09s (NLO, LO, error sets) and EKS98s (LO) are released. As an application, we consider the nuclear modification 
factor R^^f^^ at midrapidity for neutral pion production in deuteron-gold collisions at RHIC in NLO. Comparison with 
the PHENIX data in different centrality classes is also shown. In addition, predictions for the corresponding nuclear 
modification factor R"^^,^ in proton-lead collisions at the LHC are discussed. 
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1. Introduction 

In a high-energy hadronic or nuclear collision of particles A and B the inclusive cross sections for hard processes, 
where the interaction scale is large, » ^qcd' computed using the QCD collinear factorization theorem 

mm, 

^B^k^X ^ 2 jrA^Ql^ ^ j:B^q2^ ^ do-''^'^*+^' + 0(l/Q^), (1) 

where d& are the perturb atively computable partorric pieces (cross sections in lowest order), and f^^ (f^) is the parton 
distribution function (PDF) for a parton flavor / (j) in A (B). 

As is well known, the PDFs of nucleons bound to a nucleus, the nPDFs, are modified relative to the free-nucleon 
PDFs. Global DGLAP analyses, analogous to the free-proton case, have been developed also for the nPDFs. In this 
study we have considered two globally fitted sets, the EKS98 fT] and EPS09 Q, where the nPDFs are defined in terms 
of the nuclear modification Rf(x, Q^) and the free nucleon PDF f^{x, Q^) as 

ft{x,Q')^Rf{x,Q^)ff{x,Q^) (2) 



Email address: ilkka.heleniusOjyu.fi (I. Helenius) 



/Nuclear Physics A 00 (2012) i-g] 



2 



for each parton flavor /. Using the error sets of the EPS09 analysis one can also quantify the propagation of the 
nPDF uncertainties into the hard-process cross sections. So far all these globally analyzed nPDFs have been spatially 
independent. Thus, it has not been possible to compute the nuclear hard-process cross-sections in different centrality 
classes consistently with the global analyses. This problem is adressed in (I], where the spatial dependence of the 
EPS09 and EKS98 nPDFs has been determined using the framework discussed in the next section. 



2. Analysis Framework 



We introduce a nuclear modification r^{x, Q^, s) which, in addition to x and Q^, depends also on the transverse 
position s of the nucleon inside the nucleus. We define this quantity so that its spacial average yields the original 
(impact parameter independent) nuclear modification: 



Rf(x, Q^^ = jj d^s TAs)>f(x, Q\ s). 



(3) 



where Rf(x, Q^) is obtained from the global analysis EKS98 or EPS09. The key assumption is that the spatial de- 
pendence of rf{x, Q^, s) is a function of the nuclear thickness ^^(s). The functional form we adopt and test here is a 
simple power series of the thickness function. 



(4) 



The A-independent parameters d-ix, Q^) are obtained by fitting the R^{x, Q^) in Eq. jsjto the A-dependence of the 
spatially independent nPDFs. In practice, we minimize the;\f^ defined as 



Rf(x, Q^)-^JdhTA(s)rf(x,Q'-,s) 



Wf(x, 22) 



(5) 



where the weights Wf(x, Q^) are set by hand (see As can be seen from Fig.[T] we can reproduce the A-dependence 
very well with the power series ansatz when we take into account terms up to the fourth order This holds also for other 
parton flavors in the whole kinematic region considered and also for the error sets in EPS09 LO and NLO analyses. 
As an outcome of this fitting procedure we obtain the spatially dependent nPDF sets, which we refer to as EPS09s and 
EKS98s ("s" for "spatial"). These sets are now available at our webpage^ As an example, the gluon modification 
for the lead nucleus (A = 208) is shown in Fig. |2]as a function of x and s. We notice that the nuclear modifications 
at small s are larger than the spatially averaged nuclear modifications, and at larger s values (s » Ra), the nuclear 
effects die out. 



3. Applications 

As we are now equipped with these spatially dependent nPDFs, we are able to calculate the cross sections of hard 
processes also in different centrality bins. The nuclear modification factor for the production of a particle ^ in a hard 
process for a given centrality class can be calculated from (see [1 1 for details) 



R'AB(PT,yM,b2) = (6) 



where TabW is the standard nuclear overlap function, and the impact_narameter values bi and b2 for each centrality 
class can be calculated using the optical Glauber model |6|. In Fig. 3 we plot the nuclear modification factor 
for neutral pion production in d+Au collisions at RHIC at midrapidity in four different centrality bins, 0-20 %, 20- 
40 %, 40-60 % and 60-88 %. For the thickness function of deuterium, see Ref. [l]. The different lines correspond 
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Figure 1: The spatially averaged gluon modification Rg{x, Q^) at fixed 
values X = 0.001 and = 1.69 GeV^ as a function of A from the 
central sets EPS09NLO (crosses) and EPS09LO (pluses) and from the 
corresponding spatial fits EPS09sNLO (green) and EPS09sLO (blue). 
The nuclei at A < 16 (gray markers) were not used in our fits. From 1 1 1. 



Figure 2: The spatially dependent modification of gluon distribution in a 
Pb nucleus, r^^(x, Q^, s), from EPS09sNLO as a function of x and s = |s| 
at the initial scale = 1.69 GeV^. From (T|. 
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Figure 3: The nuclear modification factor /^^^^(Pr) in d+Au collisions for yijvw = 200 GeV at y = for difterent centrality classes. Calculations 
are in NLO pQCD using EPS09s and three dilferent fragmentation functions. The blue eiTor bands are computed with the EPS09s error sets and 
fDSS, and the data are from PHENIX |7 |. From [Ij. 



to NLO calculations (using INCNLO-packagq^ |[8|) with three different fragmentation functions and the uncertainty 
band is computed using the error sets in EPS09s. For the free proton PDFs we use the CTEQ6M set |9| and all the 
scales have been fixed to pj. The calculations are compared with the PHENIX data |7|, which we have multiplied 
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by a different overall factor in each centraUty bin. These factors are well consistent with the overall normalization 
uncertainties quoted by the experiment. In all centrality classes the calculations seem to agree with the data, when all 
the uncertainties in the data and in the calculation are taken into account. Especially the evolution of the slope at the 
data in the small-pj- region is well reproduced from central to peripheral collisions. 

We have also performed the corresponding calculations for the forthcoming p+Pb collisions at the LHC and these 
results are shown in Fig. |4] Due to the larger y/s^N, we get more shadowing in the small pr region than for d+Au 
collisions at RHIC. Also the basic features of our spatially dependent nPDFs are nicely visible: the nuclear effects are 
strongest in central collisions and they decrease towards peripheral collisions, and the difference between the central 
and the minimum bias collisions is rather small (see Ref. 1 1]). 
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Figure 4: The nuclear modification factor R'^^^ipj) in p+Pb collisions for i/Jnn = 5.0 TeV at y = for different centrality classes. Calculations 
are in NLO pQCD using EPS09s and three dilferent fragmentation functions. The blue eiTor bands are computed with the EPS09s error sets and 
fDSS QSJ. From [T]. 
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